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ABSTRACT 

We report our attempts to locate the progenitor of the type Ic SN 2004gt 
in NGC 4038. We use high resolution HST ACS images of SN 2004gt and have 
compared them with deep pre-explosion HST WFPC2 F336W, F439W, F555W 
and F814W images. We identify the SN location on the pre-explosion frames 
with an accuracy of 5mas. We show that the progenitor is below the detection 
thresholds of all the pre-explosion images. These detection limits are used to 
place luminosity and mass limits on the progenitor. The progenitor of SN 2004gt 
seems to be restricted to a low-luminosity high-temperature star, either a single 
WC star with initial mass > 40M Q or a low mass star in a binary. The pre- 
explosion data cannot distinguish between the two scenarios. 

Subject headings: stars: evolution — sup ernovae: general — supernovae:individual(2004gt)— 
galaxies:individual(NGC 4038) 



1. INTRODUCTION 

Stellar-evolution models predict that all stars with M Z ams ^ 8M should end their lives 
in a core-collapse induced supernova (CCSN) explosion (e.g. Heger et al. (2003); Eldridge & 
Tout (2004)). The identification of the progenitor of a CCSN just prior to explosion provides 
a critical test for models of stellar evolution. The task of finding the progenitors is compli- 
cated by the unpredictable nature of SNe. The two nearby supernovae SN 1987A (Walborn 
et al. 1989) and SN 1993J (Aldering et al. 1994) had progenitors identified shortly after dis- 
covery. Modern telescope archives have allowed the search for progenitors to be increased to 
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large distances (Maund & Smartt 2005; Van Dyk et al. 2003). In order to confidently identify 
the correct star as a SN progenitor on images of the highest resolution, precise astrometric 
techniques are required. Recently, Maund & Smartt (2005) and Smartt et al. (2004) have 
used high resolution post-explosion images with the Hubble Space Telescope (HST) and the 
technique of differential astrometry to accurately identify progenitors on pre-explosion HST 
images. The brightness and colour information from pre-explosion images may then be used 
to place the progenitor on the Hertzsprung-Russell (HR) diagram (e.g SN 2003gd, Smartt 
et al. 2004). There are no detections yet of progenitors of the hydrogen deficient SNe lb or 
Ic, although upper limits to their absolute magnitudes have been determined (see §4). It 
is now thought that some SNe Ic produce gamma-ray bursts (GRBs; Matheson et al. 2003; 
Hjorth et al. 2003), and that the progenitors may be rapidly rotating massive Wolf-Rayet 
(WR) stars. 

This letter presents the study of HST Wide Field Planetary Camera 2 (WFPC2) images of 
the site of SN 2004gt in the interacting galaxy NGC 4038, prior to explosion. SN 2004gt 
was discovered by L.A.G. Monard (Monard et al. 2004) on 2004 Dec 12.076 with an unfil- 
tered magnitude of 14.6, corresponding to an absolute magnitude of ~ —17 at the distance 
of NGC 4038 ((m - M)„ = 31.4; Whitmore et al. 1999). The position of SN 2004gt was 
given as a 2 ooo = 12 ft 01 m 50 s .37, 5 2 ooo = -18°52'12."7. Ganeshalingam et al. (2004) classified 
SN 2004gt of Type Ic. The WFPC2 data are the deepest, and widest wavelength coverage, 
pre-explosion images of any nearby SNe Ibc to date, and provide a unique opportunity to 
study the progenitor site of a possible GRB-related supernova event. Gal- Yam et al. (2005) 
use post-explosion ground-based imaging with Adaptive Optics, and similar pre-explosion 
HST WFPC2 data, and identify the same location for the progenitor of SN 2004gt and draw 
similar conclusions about its nature to those presented here. 



2. OBSERVATIONS AND DATA ANALYSIS 

The pre-explosion HST WFPC2 observations of the site of SN 2004gt were obtained 
on 1996 Jan 20 with F336W (4500s), F439W (4000s), F555W (4400s) and F814W (2000s) 
filters. The pre-explosion WFPC2 observations (Programme GO-5962; PI: B. Whitmore), in 
each filter, are drizzled combinations of four separate exposures. These four exposures were 
acquired in pairs, for the rejection of cosmic rays, with an offset of 0.25" between pairs of 
exposures (Whitmore et al. 1999). The progenitor of SN 2004gt fell on the Planetary Camera 
(PC) chip for which the pixel size of the drizzled image was 0.025" (the same as the HST 
Advanced Camera for Surveys (ACS) High Resolution Channel (HRC)). Aperture photom- 
etry was performed on the final images using the Iraf implementation of DAOphot and 
the prescription of Whitmore et al. (1999). Empirical aperture corrections were calculated 
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for each frame and small charge transfer inefficiency corrections (~ O.Olmag) were applied to 
the photometry using the prescription of Dolphin (2000a). The photometry was converted 
to the standard Johnson-Cousins system using the transformations of Holtzman et al. (1995) 
with the updates of Dolphin (2000b). 

DAOphot provides four centering algorithms to calculate the locations of stars: centroid, 
ofilter, Gaussian and PSF fitting. The scatter in the measured locations of stars was used 
to describe the uncertainty due to the under-sampled PSF of WFPC2. Post-explosion 
ACS/HRC images (Programme GO-10187; PI: S. Smartt) were acquired on 2005 May 16 
with the F435W (1672s), F555W (1530s) and F814w (1860s) filters. These images were anal- 
ysed using the PyRaf implementation of DAOphot, using Tiny Tim PSFs (Krist & Hook 
2004) for PSF photometry. Empirical aperture corrections to 0.5" were calculated for each 
frame (with the correction for 0.5" to oo taken from Sirianni et al. (2005)), and the photom- 
etry was corrected for charge transfer inefficiency according to the relationships of Riess & 
Mack (2004). The ACS photometry was converted to standard Johnson- Cousins magnitudes 
using the transformations of Sirianni et al. (2005). The photometry was found to agree with 
results from the DOLPhot package, an update of HSTPhot (Dolphin 2000b) to & 0.1 
mag atmsi 25. The locations of 16 stars common to both the pre-explosion WFPC2 imag- 
ing and the post-explosion ACS imaging were used to calculate the transformation between 
the two sets of F555W images, using the Iraf task geomap. The positions of stars on the 
ACS frame were transformed to the coordinates of the drizzled PC F555W image with an 
accuracy of ±5mas. The reverse transformation has an uncertainty of ±4.6mas. SN 2004gt 
is located 1.3", approximately South, from knot S(Whitmore et al. 1999). 



3. OBSERVATIONAL RESULTS 

SN 2004gt is readily identifiable on post-explosion ACS images. We measure the ap- 
parent brightness and colours of SN 2004gt on 2005 May 16 as my = 18.55 ± 0.12, B — V = 
0.87 ± 0.16 and V - I = 1.32 ± 0.15. The location of SN 2004gt and its position on the 
WFPC2 pre-explosion images are shown as Fig. 1. The limiting magnitude at the position 
of SN 2004gt on the pre-explosion frames was determined using the technique of Maund 
& Smartt (2005). Aperture photometry was done over a 9-point grid centered at the SN 
position, with offsets of 1 pixel which permitted background variations, over the aperture 
area, to be taken into account. The detection thresholds determined in this manner were 
tested by adding a scaled PSF at the SN position on the pre-explosion magnitude. This 
tested if the DAOphot photometry would recover a star at the calculated detection thresh- 
old. The detection thresholds were determined to be: rriFzmw ~ 23.04, mpAmw ~ 24.56, 
m F555W ~ 25.86 and m F8 i4 W 24.43. 
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The post-explosion ACS photometry of stars within 2" was used to estimate the reddening 
towards SN 2004gt using the technique of Maund & Smartt (2005). The photometry of 47 
nearby stars was used to estimate the reddening as E(B — V) = 0.07 ± 0.01. This low red- 
dening is consistent with the foreground reddening as quoted by NED 4 , after Schlegel et al. 
(1998), and the reddening determined by Whitmore et al. (1999) towards knot S. The shape 
of the continuum of an early spectrum of SN 2004gt (Kinugasa, Private Communication) 
did not show evidence for significant reddening, being similar in shape to early time spectra 
of SN 2002ap (E(B - V) = 0.09,Mazzali et al. 2002). 

The detection thresholds, calculated above, were converted to absolute magnitudes taking 
into account the distance and extinction towards NGC 4038. We calculated the extinc- 
tions in WFPC2 bands, using the A X /E(B - V) relations of Van Dyk et al. (1999), to be: 
A F33 qw — 0.39, A F439W = 0.30, A F555W = 0.20 and A F814W = 0.13. This implies the pro- 
genitor on the pre-explosion frames had magnitudes: M F336W ^ —8.75, M F439W ^ —7.14, 
M F555W £ -5.74 and M F814W £ -7.10. 

An age for the environment around SN 2004gt was estimated by comparing the positions of 
nearby stars on a colour-magnitude diagram with theoretical isochrones (we use the stellar 
evolution tracks and isochrones of the Geneva stellar evolution group 5 ). The isochrones were 
shifted for the distance, extinction and reddening of NGC 4038. Assuming the standard 
Galactic extinction-reddening relationship of Ay as 3AE(B — V) and E(V — I) = 1.6E(B — 
V), the value of the reddening calculated above yielded Ay = 0.22 ± 0.03 and E(V — I) = 
0.11±0.05. The mean age of the nearby stars was measured to be log(age/years) = 6.93 ± 0.13, 
or 8.5 ± 2.5Myr. This age is consistent with the age of 7 ± lMyr estimated by Whitmore 
et al. (1999) for knot S. The age of the nearby stars corresponds to the lifetime of stars with 
M ZA ms ~ 20 - 4OM . 



4. DISCUSSION 

Maund & Smartt (2005), Smartt et al. (2002) and Van Dyk et al. (2003) have shown 
how the pre-explosion photometry may, in the event of the non-detection of the progenitor, 
be used to place mass limits on the progenitor. The pre-explosion detection thresholds were 
converted to luminosity thresholds for a range of supergiant spectral types, using the colours, 
temperatures and bolometric corrections given by Drilling & Landolt (2000). In addition 
synthetic photometry was conducted on model WR spectra (Grafener et al. 2002; Maund & 
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Smartt 2005), using the STSDAS packages synphot, to calculate the colours and bolometric 
corrections for these types of stars. These colours and corrections are dependent on two 
parameters: the effective temperature and the "Transformed Radius." Maund & Smartt 
(2005) discuss how the two-parameter space of model spectra was sampled to completely 
measure the spread in colours and bolometric corrections. The principal consequence of this 
two-parameter dependence is that at high temperatures, when the progenitors are expected 
to be WR stars, the uncertainty in the luminosity increases with temperature compared 
to the relationship for supergiants. The observed absolute magnitude detection limits were 
converted to luminosities using Eqtns. 6 and 7 of Maund & Smartt (2005). The detection 
thresholds for the four different filters were plotted on HR diagrams, shown in Fig. 2, and 
compared with stellar evolution tracks to estimate mass limits. 

The F336W pre-explosion observation does not place a stringent constraint of the progenitor, 
missing both the red supergiant and blue WR end points of the stellar evolution tracks 
for low-intermediate and high-mass stars, respectively. The F439W observation excludes 
yellow and red supergiant progenitors, for Mzams ~ 20 — 25M Q as well as the lowest mass 
stars which end their lives with a WR phase (M ZAMS 25 — 4OM ). The constraint 
against a WR progenitor, with initial mass of 4OM , is weak , however, as the luminosity 
uncertainty is large. The red supergiant progenitors for stars with Mzams < 25M are 
completely excluded by the F814W observation (including the location of the progenitor of 
the type IIP SN 2003gd; Smartt et al. 2004). The F555W observation places much stronger 
constraints on both sides of the HR diagram. The red supergiant progenitors for stars with 
Mzams < 25M Q are completely excluded (including the location of the progenitor of the type 
IIP SN 2003gd; Smartt et al. 2004). A much tighter constraint on the WR branch excludes 
a Mzams ~ 25 — 4OM progenitor. All four pre-explosion observations clearly indicate that 
the SN 2004gt was not the explosion of a super-massive M Z ams ~ 12OM star. The depth 
of the observations, despite not being sufficient to detect the progenitor, provide useful limits 
on the progenitor of SN 2004gt. The use of the smaller distance of 13.8 ± 1.7Mpc to NGC 
4038, determined by Saviane et al. (2004), would result in the threshold limits on the HR 
diagrams being lowered by ~ O.281ogL . The consequence of this would be to increase the 
likelihood of the detection of a WR star progenitor arising from stars with initial masses 
~ 4OM . This shift is, however, smaller than the uncertainties of the thresholds at the blue 
side of the HR diagram. The pre-explosion observations do not exclude possible blue lower 
mass progenitors in binaries, with Mzams = 20 — 4OM as given by the age determined 
for the surrounding stellar population. The evolution of such a progenitor would have been 
affected by angular-momentum and mass exchange with a companion, which would strip 
the progenitor's H envelope leaving a low-luminosity high-temperature He and C-0 star 
as the progenitor (Podsiadlowski et al. 2004). Maund & Smartt (2005) have suggested 
this as a possible explanation for the lack of the detection of a number of progenitors of 
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type Ib/c SNe. This progenitor scenario was also invoked by Nomoto et al. (1994) for the 
type Ic SN 19941. Podsiadlowski et al. (2004) argue that prompt collapse of a single star 
would produce a faint SN, whereas SN 2004gt is of normal brightness. The distance to 
NGC 4038, 19.2Mpc, is prohibitive to conducting a spectroscopic search for the companion 
(Maund et al. 2004). In Table 1 we compile all the directly measured upper limits for the 
magnitudes of the progenitors of SNe Ibc, in various pass-bands. We can determine the upper 
final mass limit, i.e. the mass prior to explosion, of a WR star progenitor from the mass- 
luminosity relationship of the Langer models(Langer 1989). These upper limits are somewhat 
higher, although not inconsistent with the calculated final mass of the star which produced 
SN 2002ap (5 — 7M Q ; which includes ~2M allowance for a compact remnant) determined 
from modelling the spectral evolution of the supernova (Mazzali et al. 2002). Vacca & 
Torres-Dodgen (1990) give the range of observed absolute magnitudes for Galactic WN stars 
as My = —3.2 to 8.0. Assuming a uniform distribution in absolute brightness implies that 
the pre-explosion observations are sensitive to 50% of WN stars, mostly WNL stars. 
The luminosities of WC stars in the LMC (Crowther et al. 2002) are log(L/L ) < 5.8, which 
places them within the lower uncertainty boundary for the detection of WR stars (see Fig. 2). 
This shows that these observations cannot constrain the progenitor scenario by themselves. 
A comparison of these final upper mass limits may, however, be used with SN evolution 
models and their mass budget to constrain a WR progenitor scenario. Understanding the 
progenitors of Type Ic SNe, whether high-mass WR stars or low-mass stars in binaries, is 
important for understanding the relationship between SNe and GRBs (Hjorth et al. 2003). 
The high rate of SNe events in the merging galaxies NGC 4038/4039 (0.2 — 0.3 SNe yr _1 ; 
Neff & Ulvestad 2000) in conjunction with the large amount of deep imaging of this system 
available suggests that NGC 4038/4039 is a good candidate for the detection of progenitors 
of Type Ibc SNe in the future. 

The location of SN 2004gt, in the galaxy NGC 4038, in pre-explosion HST WFPC2 
images has been identified using differential astrometry with a post-explosion high-resolution 
HST ACS observation of SN 2004gt. The progenitor location was identified with an accuracy 
of 5mas. The progenitor was not detected in any of the deep F336W, F439W, F555W and 
F814W pre-explosion imaging. The detection thresholds of the pre-explosion observations 
suggest that SN 2004gt arose from a single star progenitor with M Z ams ^ 40M Q or a 
Mzams ~ 20 — 4OM star in a binary. This is the most restrictive mass limit placed on 
any type Ib/c SN yet, and this is attributable to the depth of the available pre-explosion 
imaging. 

Based on observations made with the NASA/ESA Hubble Space Telescope and obtained 
from the Data Archive at the Space Telescope Science Institute, which is operated by the 
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Association of Universities for Research in Astronomy, Inc., under NASA contract NAS 5- 
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NSF through grant AST02-05994. 
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Table 1: Compilation of all the directly determined limits to the progenitors of SNe Ibc. 



Supernova 


Type 




Limiting 


; mag 
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log L ho i 
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M B 


My 


Mr 


Mi 


(L Q ) 


(M ) 


2000ds 


lb 






-6.0 




-6.5 


5.2 


10 


2000ew 


Ic 






-6.4 






5.5 


15 


2001B 


lb 






-8.0 






6.0 


25 


2002ap 


Ic 


-8.7 


-6.8 


-7.2 


-7.8 


-8.4 


5.6 


18 


2004gt 


Ic 


-8.42 


-6.82 


-5.31 




-6.77 


5.1 


9 



NOTES. — Limits are from this paper, Smartt et al. (2002), and Maund & Smartt (2005), 
with the corresponding bolometric luminosity limit calculated from Maund & Smartt (2005). 
The limiting magnitudes for SN 2002ap have been adjusted for the new large distance to 
M74 as determined by Hendry et al. (2005). The M v and Mj magnitudes for SN 2000ds 
are HST M F555W and M F8 uw] My for SN 2000ew and SN 200 IB are M Fm6W and M F555W 
respectively. For SN 2004gt the limiting magnitudes M v , M B , My and Mj are HST M F336W , 
M F439W , M F555W and M F8UW respectively. 
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Fig. 1. — Pre- and Post-explosion images of the site of SN 2004gt in NGC 4038. The same 
orientation is used for all three panels, with North up and East to the left. SN 2004gt is 
located at the centre of each of the panels, a) Post-explosion ACS/HRC F555W image, with 
SN 2004gt clearly visible, b) Post-explosion ACS/HRC F814W image, c) Pre-explosion 
WFPC2 PC chip F555W image, transformed, shifted and rotated to the coordinates of the 
post-explosion ACS imaging. No object is significantly detected at the position of SN 2004gt 
prior to explosion, with a positional uncertainty 0.004". In these images knot S appears in 
the top-right. 
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Fig. 2. — Hertzsprung- Russell diagrams showing solar metallicity stellar evolution tracks. 
Overlaid are the detection thresholds for the pre-explosion HST WFPC2 F336W, F439W, 
F555W and F814W images. A progenitor lying in the grey regions would have been signifi- 
cantly detected. 



